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Abstract 

In 1,2 or 3 dimensions a scalar wave excited by a non-negative source 
in a viscoelastic medium with a non-negative relaxation spectrum or a 
Newtonian response or both combined inherits the sign of the source. 
The key assumption is a constitutive relation which involves the sum 
of a Newtonian viscosity term and a memory term with a completely 
monotone relaxation kernel. In higher-dimensional spaces this result holds 
for sufficiently regular sources. Two positivity results for vector-valued 
wave fields including isotropic viscoelasticity are also obtained. 



Notation. 

[a,b[:= {x eR \ a < X <b}; 
R+ = ]0,cx)[ ; 
I: unit matrix; 
(k,x) := J2t=i hx^] 

/,/ N f 1, y > 



Hp) :=/(re-P^/(y)dy. 
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1 Introduction 



Positivity of viscoelastic pulses was studied in a paper oflP^ (|l969l) . Duff 
assumed a special model with a rational complex modulus. Duff's models are 
however loosely related to viscoelasticity and his assumptions are excessively 
restrictive. 

In this paper a general scalar viscoelastic medium with the constitutive equa- 
tion a = a e + G(t) * e with a completely monotone relaxation modulus G and a 
non-negative Newtonian viscosity coefficient is studied. We show that a scalar 
viscoelastic wave field propagating in a d-dimensional medium and excited by 
a non-negative pulse is also non-negative provided d < 3. For higher dimen- 
sions and for non-zero initial data only wave fields excited by sufficiently regular 
sources are non-negative. 

Positivity of viscoelastic signals can be considered as a test for the non- 
negative relaxation spectrum and for the presence of the Newtonian viscosity. 

Positivity can be extended to matrix-valued fields, e.g. to Green's functions 
of systems of PDEs. In Sec. [H] we consider a system of PDEs resembling the 
equations of motion of viscoelasticity with a CM relaxation kernel and prove that 
the Green's function of this system of equations is positive-semidefinite. This 
result does not apply to general viscoelastic Green's functions, which involve 
double gradients of positive semi-definite functions. In iso 



2 Statement of the problem 

In a hereditary or Newtonian linear viscoelastic medium a scalar field excited 
by positive source is non-negative. This applies to displacements in pure shear 
or to scalar displacement potentials. The key assumption about the material 
properties of the medium is a positive relaxation spectrum. The result holds for 
arbitrary spatial dimension. 
We consider the problem: 

pD^u^aV^Du + G{t)*W'^'Du + s{t,x) t>0, x eR'^ (1) 

with s(i, x) — 9{t) (ci + C2 t) 6{x) and the initial condition 

u{0,x) — uq S{x), 'Du{0, x) = uq 6{x), (2) 

(Problem I) as well as s{t, x) — c5{t) 5{x) with a solution assumed to vanish for 
t < (Problem II). It is assumed that a > and G is a completely monotone 
(CM) function. 

The Laplace transform 

POO 

u{p,x):^ e-P^u{t,x)dt, Rep > 0, x eR'^ (3) 

satisfies the equation 

pp^ u(p, x) = Q{p) u(p, x) + g{p) S{x) (4) 



2 



where 

Q{p):^ap + pGip) (5) 
The function g is defined by the equation 

9i{p) =^ ^+ puo + iiQ (6) 

in Problem I and 

91 lip) = 1 (7) 

in Problem II. 



3 Basic mathematical tools. 

The clas ses of functions appropriate for viscoelastic responses are reviewed in 
detail in Seredvnska fc Hanygal ( 20091 ). 

Theorem 3.1 // the function u{-, x) is completely monotone for every x G , 
then u{t, x) > for every t > and x £ R*^. 

Proof. li u(-,x ) is completely monotone, then, in view of Bernstein's theorem 
( Widder . for every x S R'' it is the Laplace transform of a positive Radon 

measure rrix- 



u{;x) = / e-P'm^ds) (8) 

The Radon measure is uniquely determined by u{-^x), hence mx{dt) = 
u{t,x) dt is a positive Radon measure. Hence, in view of continuity of u{-,x), 
we have the inequality u{t, x) > for alH > and x £ R''. □ 

The problem of proving that u{t, x) is non-negative is thus reduced to proving 
that u{-,x) is completely monotone. The crucial step here is the realization 
that Q in ([5]) is a complete Bernstein function. We shall therefore recall some 
facts about Bernstein and complete Bernstein functions and their relations to 
completely monotone functions. 

Definition 3.2 A function f on R+ is said to be completely monotone (CM) 
if it is infinitely difjerentiable and satisfies the infinite set of inequalities: 

(—1)" D" /(y) > y > 0, for all non-negative integer n 

It follows from the definition and the Leibniz formula that the product of two 
CM functions is CM. A CM function can have a singularity at 0. 

Definition 3.3 A function f on R+ is said to be locally integrable completely 
monotone (LICM) if it is CM and integrable over the segment ]0, 1]. 
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Definition 3.4 A function f on R+ is said to be a Bernstein function (BF) if 
it is non-negative, dijferentiable and its derivative is a CM function. 

Since a BF is non-negative and non-decreasing, it lias a finite limit at 0. It can 
therefore be extended to a function on M+. 

Every CM function / is the Laplace transform of a positive Radon measure: 



Theorem 3.5 (Bernstein's theorem Widdei i '19AA )) 



10,' 



(9) 



It is easy to show that / is a LICM if the Radon measure n satisfies the inequality 

fi{dr) 



[0,' 



1 + r 



< oo 



(10) 



Theorem 3.6 UacoA . l20M : \Seredvnska & Hanva2. \200A ) If f,g are CM then 
the pointwise product f g is CM. 

Let / be a Bernstein function. Since the derivative D/ of / is LICM, Bern- 
stein's theorem can be applied. Upon integration the following integral repre- 
sentation of of a general Bernstein function / is obtained: 



f{y)=a + by+ f [l-c-'^y] iy{dr) 

JlO,oo[ 



(11) 



where a,b = D/(0) > 0, and i^(dr) :— fi{dr)/r is a positive Radon measure on 
M+ satisfying the inequality 



0,oo[ 



r iy{dr) 

1 + r 



< oo 



(12) 



The constants a, b and the Radon measure v are uniquely determined by the 
function /. 

Theorem 3.7 iBera & Forsi \l97d : \Jacol]i. \200i ) If f is a CM function, g is 
a BF and g{y) > for y > then the composition f o g is a CM. 

Corollary 3.8 \Berg & Forsi \l97d : \jaco^ . \200i ) If g is a non-zero BF then 
1/g is a CM function. 

Note that the function f{y) := exp(— y) is CM but 1// is not a BF. 

Definition 3.9 A function f is said to be a complete Bernstein function ( CBF) 
if there is a Bernstein function g such that f{y) — y^ g{y). 
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Theorem 3.10 l(Jacol\ . \200l ) A function f is a CBF if and only if it satisfies 
the following two conditions: 

1. f admits an analytic continuation f{z) to the upper complex half-plane; 
f{z) is holomorphic and satisfies the inequality Im/(z) > for hnz > 0/ 

2. f{y)>OforyeR+. 

The derivative Dg of the Bernstein function 5 is a LICM function h. Hence 
we have the following theorem: 

Theorem 3.11 Every CBF f can be expressed in the form 

f{y) = yh{y) + ay (13) 

where h is LICM and a = g{Q) > 0. Conversely, for every LICM function h 
and a > the function f given by (jlSp is a CBF. 



Proof. For the first part, let g be the BF in Definition 13 . 91 and let h := Dg. Since 
h{x) dx = g{l) — g{0) < 00, the function h is LICM. For the second part, 
note that if h is LICM, then g{y) — a + h{s) ds is a BF and f{y) = y^ g{y)- 

a 

Since the Laplace transform of a LICM function h has the form 

Hv) = f ^ (14) 

J[o,oc.[ r + y 

where fj, is the Radon measure associated with h, every CBF function / has the 
following integral representation 

f{y)^b + ay + yf ^ (15) 

^]o,oo[ r + y 

with arbitrary a,b ~ m({0}) > and an arbitrary positive Radon measure fj, 
satisfying eq. (|10p . The constants a,b and the Radon measure /i are uniquely 
determined by the function /. 

Noting that y/{y + r) = r [1/r — 1 / {y + r)], we can also express the CBF / 
in the following form 



/•oo 

f{y)^b + ay+ [l-e-'y]h{z)dz 
Jo 

where h{z) :— J^^ e"*"^ m(dr) > and m(dr) := r/i(dr) satisfies the inequal- 
ity 

f m(dr) 

/ , ' < 00 

J[o.oo[ r(r + 1) 
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Let i'{dz) h{z) dz. Wc have 



zv{dz) f°°zh{z)dz 



[0,oo[ 1 + 2: Jo l + Z J[Q 



m{dr) [l/r-e''r(0,r)] 



,oc 



Using the asymptotic properties of the incomplete Gamma function (jAbramowitz fc Stegun 



19701) it is possible to prove that the right-hand side is finite, hence the Radon 



measure i/{dz) := h{z)dz satisfies inequahty (fT2| . We have thus proved an 
important theorem: 

Theorem 3.12 Every CBF is a BF. 

However 1 — exp(— y) is a BF but not a CBF. 
The simplest example of a CBF is 

Va{y) ■■= yl{y + a)=y' e"^^ [1 - e-^''] ds 
Jo 

a > 0. It follows from eq. that every CBF / which satisfies the conditions 
/(O) = and lim^^oo fiy)/y = is an integral superposition of the functions 
(pa- The CBF ipa corresponds to a Debye element defined by the relaxation 
function Ga{t) = exp(— a<). 

We shall need the following properties of CBFs: 



Theorem 3.13 UacoA. \200l : \Seredvnska & Hanvaci . 1200. 



1. f is a CBF if and only if y/ f{y) is a CBF; 

2. if f, g are CBFs, then f o g is a CBF. 

The second statement follows easily from Theorcm l3.10l 

Remark. 3.14 is a CBF i/0 < a < 1, because 

1 r°° 

?y""^ = / e~^*s~"ds = 

^ r(i-a)yo 



00 /'OO 



— \ , , / dz I e-y' e-^'ds z"-^ 

r(l - a) r(a) Jo Jo 



' ■ dz 



7t Jo y + z 
and thus 

^ y sin(a7t) 

Jo y + z 

The sets of LICM functions and CBFs will be denoted by ^ and €. respec- 
tively. 
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4 Positivity of one- and three-dimensional solu- 
tions. 



Applying the results of the previous section, we get the following result: 

Theorem 4.1 If a > and the relaxation modulus G is CM then the Junction 
Q defined by eq. ^ is a CBF. 

The mapping {a,G) G x ^ Q E € defined by eq. ([5]) is bijective. 

A one-dimensional solution of eq.([4]) is given by 

ui{p,x) = Ui{p, \x\) := A{p) exp(-B(p) |x|) 

with B{p) = pi/2p/Q(p)i/2 and A{p) = g{p)/[2B{p)] If Q g £, then Q{yY'^ 
is a composition of two CBFs, namely y^^^ fRemark I3.14p and Q, hence it is 
a CBF by Theorem EII3 The function B{p) is a CBF by Theorem OS] and 
l/B{p) is a CM function by Theorem [SH and Corollary EjH 

The amplitude of the solution of Problem I is given by A{p) = l/\2p B{p)] + 
uo/[2i?(p)]. The first term is a CM function because it is the product of two 
CM functions. The second term is also CM, hence A{p) is CM. The amplitude 
of the solution of Problem II A{p) = l/[2 B{p)\ is also CM. 

For every fixed x the function exp(— i3(p) \x\) is the composition of a CBF 
and the function B, which is a CBF and therefore a BF. By Theorem 13.71 the 
function exp(— |a;|) is CM. This proves that for d = 1 the solutions of 
Problem II and Problem I with mq = arc non-negative. 

In a three-dimensional space the solution of ^ is given by the equation 

1 dU^jp^r) 
-3(p,a.) = - — 

where r — \x\, so that 

U3(p, x) = -^A{p) B{p) exp{-B{p) \x\) (16) 

But A{p) B{p) — g{p)/2. If mq = then g is CM. Hence u^{-,x) is the product 
of two CM functions and thus CM. 



5 Positivity of solutions in arbitrary dimension. 

In an arbitrary dimension d 



Ud{p,x) 



9{p) 



{2nYQ{p) 



1 



(17) 



The above formula ca n be expressed in terms o f MacDonald functions by using 
eq. (3) in Sec. 3.2.8 of ICeFfand &: Shilovl (|l964 : 



Ud{p,x) 



P' 



,d/4-l/2 



9ip)P 



,d/2-l 



(27t)''/2 Q(p)'i/4+l/ 



— r-W2-i)if,/2_i(B(p)r) 



(18) 
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where B{p) is defined in the preceding section. 

The MacDonald function is given by the integral representation 

/•OO 

Kfj^{z)= I exp(— z cosh(s)) cosh(/is) ds (19) 
Jo 

Since cosh(?/) is a positive increasing function, it follows immediately that Kf^{z) 
is a CM function. 

We shall need a stronger theorem on complete monotonicity of MacDonald 
functions. 

Theorem 5.1 \Miller & Samk^ . \200i ). 

The function z^/^ K^{z) is CM for > 1/2. Q 

The proof of this theorem requires a lemma. 

Lemma 5.2 If a > then the function (1 + l/x)" is CM. 
Proof. We begin with < a < 1. Setting t — l/{xy) we have that 

' ^ " ' dt = 



Ji (xy + 1)1-" Jo (l/i+l)!-" 

a J + - 

Since for each fixed value oi y > the function {xy + 1)"^^ is CM, the function 
(1 + l/x)" {x > 0) is also CM. 

The function 1 + 1/a; is CM, hence for every positive integer n the function 
(1 + 1/a;)" is CM. We can now decompose any positive non-integer a into the 
sum a ^ n + (3, where n is a positive integer and < (3 < 1. Consequently 

(1 + 1/x)" = (1 + 1/x)" (1 + 1/xf 

is CM because it is a product of two CM functions. □ 

Proof of the theorem For /i > —1/2 the MacDonald function has the following 
integral representation: 

' V 2r(i/2-A*) 7o ^ 2zy 

, , , ^ (20) 

(jCradshtevn fc RvzhikI (|l994[ ). 8:432:8). By Lemma [Ol the integrand of the 
integral on the right-hand side is CM ii ^ > 1/2. Hence the integral is the 
limit of sums of CM functions, therefore itself a CM function. Consequently, 
the function z^/"^ Kfj,{z) is the product of two CM functions, and thus it is CM 



The theorem is vaUd for /i > 0, see iMiller fc Samkol ll200lf) . but we do not need this fact. 
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too. 



□ 



We now note that u{p,x) = p^'^~^^^^ g{p) F{p). We shall prove that F{p) 
is the product of two CM functions of the argument p, viz. Q{p)^'--'^^^^^'^ and 
L{z) := z^/^ ^(j/2-i(^) with z := B{p) r, as well as a positive factor independent 
of p. 

Lemma 5.3 If Q is a CBF and a > 0, then Qipy" is CM. 

Proof. Let n be the integer part of a, a = n + /?, > /3 < 1. Q{p)~^ is 
CM (by Theorem and Corollary [Sj]) and therefore also is CM. By 

Theorem I3TT31 the function Q{p)'^ is a CBF, hence l/Q{p)'^ is CM. Consequently 
Qip)-" is CM. □ 
The lemma implies that the factor is CM. Since the function L 

is CM and we have already proved that B{p) is BF, Theorem 13.71 implies that 
L{B{p) r) is a CM function of p. For d < 3 the factor is also CM. 

Consequently, for d < 3 the solution u(t^ x) of Problem II is non-negative. The 
solution of the same problem with an arbitrary source of the form s{t) 6(x) and 
s{t) > can be obtained by a convolution of two non- negative functions and 
therefore is also non-negative. 

For d < 5 Problem I with uq — uq — has a non-negative solution if ci > 0. 
For d <7 Problem I has a non- negative solution if ci = and C2 > 0. 

For d > 3 the fractional integral 

1 /■* 

P u{t, x)^—- {t~ f{s) ds, t>0 

Jo 

is non-negative provided a> (d — 3)/2 and uq ~ or provided a > {d — l)/2. 
We summarize these results in a theorem. 

Theorem 5.4 In a viscoelastic medium of dimension d < 3 with a constitutive 
relation 

a = ae + G(t)*e, a > 0; Ge^ 

Problem II as well as Problem I with the initial condition uq = have non- 
negative solutions. 

Under the same assumptions but for an arbitrary dimension d > 3 certain 
indefinite fractional time integrals of the solution are non-negative. For zero 
initial data Problem I has a non-negative solution if d < 5 and C\ > 0, or if 
d < 7, ci = and C2 > 0. 

6 Positivity properties of vector-valued fields. 

It is interesting to examine the implications of CM relaxation kernels on positiv- 
ity properties of vector fields. We shall prove that in a simple model complete 
monotonicity of a relaxation kernel implies that the Green's function is positive 
semi-definite. 
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Unfortunately the tools developed in Sec. [3] fail for matrix- valued CM and 
Bernstein functions Q{p) which do not commute with their derivatives. In 
particular, the product of two non-commuting matrix-valued functions need 
not be a CM function and the function / o G, where / is CM and G is a 
matrix-valued BF, need not be CM. 

Definition 6.1 A matrix-valued function F : R+ M"^" is said to be a CM 
function if it is infinitely differentiable and the matrices (—1)" D" F{y) are pos- 
itive semi-definite for all y > 0. 

Definition 6.2 A matrix-valued Radon measure M is said to be positive if the 
matrix (v, /jp^^j /(y) M(d?/) v) > for every vector v g M" and every non- 
negative function f on with compact support. 

It is convenient to eliminate matrix-valued R adon measures by applying the 



following lemma (jHanvga fc Seredvhskal . 120071 ): 



Lemma 6.3 Every matrix-valued Radon measure M has the form M(da;) = 
K(x) m{dx), where m is a positive Radon measure, while K is a matrix-valued 
function defined, bounded and positive semi-definite on except on a subset 
E such that m{E) — 0. 



Theorem 6.4 ICripenbera et al" . 199 A ) A matrix-valued function F : 



K"^" is CM if and only if it is the Laplace transform of a positive matrix-valued 
Radon measure. 

The following corollary will be applied to Green's functions: 

Corollary 6.5 //R(p) :== e-^* R(t) dt is a matrix-valued CM function then 
R(t) is positive semi-definite for t > 0. 

Definition 6.6 A matrix-valued function G : R"^" is said to be a 

Bernstein function (BE) if G{y) is differentiable and positive semi-definite for 
all y > and its derivative DG is CM. 

Definition 6.7 A matrix-valued function H : R+ R"^" is said to be a 
complete Bernstein function (CBF) if 'H.{y) = y'^G(y), where G is an n x n 
matrix-valued BF. 

The integral representation (|15p of a CBF remains valid except that the 
Radon measure has to be replaced by a positive matrix-valued Radon measure 
N(dr) = K(r)z/(dr): 

H(,).B + ,A + ,/ (21) 

J]o,oo[ r + y 
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where the Radon measure v satisfies the inequality 



/ < (22) 

the matrix- valued function K(r) is positive semi-definite and bounded i^-almost 
everywhere on IR+ while A, B are two positive semi-definite matrices. Every 
matrix-valued CBF H can be expressed in the form 

H(y)=yF(y) + yA (23) 

where F is a matrix-valued LICM function. 
We now consider the following problem 

pD^g = ADg + G* V^Dg + (5(t)(5(x)I, t > 0, x e M'' (24) 

where A is a positive semi-definite n x n matrix and G is an n x n matrix- valued 
relaxation modulus. 

If the relaxation modulus G is a CM matrix-valued function then the func- 
tion Q(p) := pG{p) is a matrix-valued CBF. The function Q is real and positive 
semi-definite, hence it is symmetric and has n eigenvalues Qnip) and n eigenvec- 
tors Bk, k = 1, . . . ,n. We shall now assume that the eigenvectors are constant: 

n 

Q{p) = ^ qkip) efc » Bfc 
fc=i 

It is easy to see that the functions Qk, k = 1, . . . ,n, are CBFs. 

The Laplace transform Q {p, x) of the Green function is given by the formula 

GiP^ ^) = / ^'^•'^^^ [P' I + |kP Q(P)] d,k EE 

n „ 

J^^W)~'J ''^'"^ ''''' ' ® " 

n 

^gk{p) efc Gfc 

fc=i 

where 

d/4-1/2 „d/2-l 

^'^^ = (2n)^/^,fc(p)^/^+V2 (^^(^) ^) 

and Bk{p) p^^'^ p/qkipY^^, k — 1, . . . ,n. Assume for definiteness that d < 3. 
The argument of Sec. [5] now leads to the conclusion that the functions gk, 
fc = 1, . . . , n, are CM, hence the function t/(-, a;) is a matrix-valued CM function 
and therefore the Green function Q{t,x) is positive semi-definite for t > 0, 
X G M.'^. In particular, we have the following theorem: 
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Theorem 6.8 Let p G R+, d<3, s{t, x) = S{t) S{x) w, where w G M". 

IfG{s) = Gk (s) e/c (8) Bk and A = (gi e/c wii/i CM functions 

Gk and real numbers Ofc > 0, fc = 1, . . . , n, iften the solution u o/ t/ie problem 

pD^u^^ ADu + G* V^Du + s(t,a;), t>0, xeR"^ 

satisfies the inequality 

{u{t, x), w> > 0, t>0, xe M*^. (25) 



7 Positivity in isotropic viscoelasticity. 

Consider now the Green's function Q oi a 3D isotropic viscoelastic medium. The 
function Q is the solution of the initial- value problem: 

P Qkr {t,x) = G{t) klmn * D Gmr.nl + S{t) S {x) 5kr , 

t>0, xeR^, k,r = 1,2,3 (26) 

with zero initial conditions, and 

Gklmn (t) = A(i) 6kl 5mn + P{t) 5km ^In + Skn ^Im (27) 

where the kernels X{t),iJ,{t) are CM and p G R+. The function G with the 
components Gkimn takes values in the linear space S of symmetric operators 
on the space S of symmetric 3x3 matrices. It is easy to see that under our 
hypotheses this function is CM: 

(-1)" (ei,G(t)e2) >0 for all n = 0,1, 2... 

for every ei, 62 G S, where (v, w) := Vki Wki is the inner product on S. 
The Laplace transform 5 of 5 is given by the formula 

g(p, x) = — {V^ VA-i Fl(p, |a;|) + [l - V O VA"!] F^ip, \x\)} 
pp 



where A := V^, 



and 



JL(p,r):=%^e-f'''«^(f)'- (28) 
47rr 

j.^(p,,):=flMle-^'^'^-W'- (29) 
4nr 
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Since qh(,p) = p/sh{p) and qT{p) = p/st{p) are CBFs, the functions p/qhip)^^"^ = 
P^^'^ si,{p) and p/qiipY^"^ = p^^"^ st{p) are BFs. hence the exponentials in 
eqs (PSH^ are CM functions of p. Moreover the functions sl(p)^ and St(p)^ 
are CM. It foUows that the functions Fi,{p,r) and Ft(p, r) are CM and there- 
fore they are Laplace transforms of non-negative functions Fi^{t^ r) and F'T{t, r). 
Their indefinite integrals fi^{t,r) := Fi^{s,r)ds and /T(i,r) := /p Fi^{s,r)ds 
are also non-negative. The functions hi^{t,r) fi,{t,r) ds, hT{t,r) = 

fTit^r) involve a convolution with a non-negative kernel and therefore are 
non-negative. The Green's function can be expressed in terms of these functions: 

g{t,x) ^ i{V®V/iL(i, |a;|) + [AI- V®V] hT{t,\x\)} (32) 
P 

We shall use the notation v > if Wfc > for k = 1,2,3. 

Theorem 7.1 Let u = V(/) + V x be the solution of the initial-value problem 

pD^ u = G *\/^Du + s{t,x), t>0, xeR"^ (33) 

with u(0, x) — — Du(0, x) and s{t, x) = V/(i, x) -|- V x g{t, x) . 

Then Vf{t, x) > Q for all t > 0, x e implies that V4>{t, x) > Q for all 
t > 0, X e 

Similarly, V x g(i, x) > for allt>0,x£ implies that V x 'i/;(t, x) > 
for allt>0,x€ R^ . 

Proof. Substitute u = V0 + V x ^, s = V/ -I- V x g in the formula 

u{t,x)=J J g{t - s,x - y)s{s,y)d:ixds 

where Q is given by p2p . Noting that A^^ is a convolution operator commuting 
with V and V div s = V/ we have 

V0(i, x) = - f fi^it -s,\x- y\) (V/)(s, y) dgy 
P J 

We now note that [I — A~^ V ® V] s = s — V/ = V x g. Hence 

V X V(i, x) = -[ frit - s,\x- y\) V x g(s, y) dsy 
P J 

The functions /l and /t are non-negative, hence the thesis follows. 

□ 



8 Concluding remarks. 

A non-negative source term excites a non-negative viscoelastic pulse. This result 
holds for scalar waves and for scalar potentials under the usual assumption that 
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the stress response is determined by a CM relaxation modulus G or by a New- 
tonian term or both connected in paraUel. The CM property of the relaxation 
modulus is a fairly general property of real viscoelastic media, equivalent to the 
assumption that the relaxation spectrum is non-negative. A generalization of 
positivity for vec tor-valued visc oelastic fields in viscoelastic media with the P 
class anisotropy ( Hanvea . 20031) is sketched. 

A particular example of a CBF is the rational function F{p) = Rn{p)/ Sm{p), 
where Rn and Sm are two polynomials with simple negative roots Xk, k = 
1, . . . , N, fii, I — 1, . . . , M, M = N or N + 1 satisfying the intertwining condi- 
tions: 

< Ai < ^1 < . . . /lAT [ < Xn+1 ] 

(the last inequality is applicable only if M = iV -I- 1) (DuS, 1969f) . A more 
general CBF is obtained by substituting in F the CBF p", with < a < 1: 



F^{p)^RNip'')/SMip") 



(Theorem 13. 13|) . The choice oi Q = Fa corresponds to a generali zed Cole-Cole 
model of relaxation. For N = M = 1 the original Cole-Cole model ( Co\e fc Cold . 
Il94ll: iBaglev fc Torvikl . Il983l ) is recovered. 

Anisotropic effects can be introduced by replacing the operator by g'^' dk di . 
Hhkig'"" = 6'^ then 



Ud{p, x) = y/detg ^ 



d/4-1/2 



gip) p' 



,d/2-l 



{2n)d/2Q{p)<i/4 



—r-^'^n-i) K,„_,[p^/^pr/Q{pY/^) 



(34) 



If Q i s a CBF then u{t,x) > 0. with r := [hu x'' x^]^^^ , cf iGel'fand fc Shilovl 
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